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MANNED MARS LANDER LAUNCH-TO-RENDEZVOUS ANALYSIS 

FOR A 1981 VENUS-SWINGBY MISSION 

By Nickolas L. Faust and Thomas B. Murtagh 
Manned Spacecraft Center 

SUMMARY 


A manned lander vehicle may someday be used to research and explore the Mars 
surface. The return of this lander by way of a time-optimum launch from the Mars 
surface to some intermediate orbit, with subsequent maneuvers (calculated by using 
Lawden's primer-vector theory in an optimum multi-impulse rendezvous program) to 
rendezvous with the primary spacecraft in its Mars parking orbit is the subject of this 
report. The type of Mars mission used to demonstrate the analytical technique in- 
cludes a Venus swingby on the Mars-to-Earth portion of the trajectory in order to re- 
duce the total mission velocity requirement. The total velocity requirement, assuming 
inplane launches, was approximately 5334 m/sec (17 500 ft/sec) for the nominal 30-day 
stay on the planet surface. 


INTRODUCTION 


One type of Mars mission includes a Venus swingby on the Mars-to-Earth por- 
tion of the trajectory in order to reduce the total mission velocity requirement. The 
stay time in Mars orbit for this mission is 80 days with 30 days nominally planned for 
research and exploration of the Mars surface. A manned Mars lander will be used to 
accomplish this major scientific objective. The return of this lander to the primary 
spacecraft by use of a launch from the surface to some intermediate orbit, with sub- 
sequent maneuvers to rendezvous with the primary spacecraft (called the orbiter) in a 
Mars parking orbit, is the subject of this report. 

The analysis is presented in three phases. The first phase consists of the cal- 
culation of time-optimum, constant-thrust, two-dimensional launch profiles to the 
intermediate orbit. Results of a parametric analysis of these launch profiles indicate 
the appropriate lander characteristics as well as the shape of the intermediate orbit. 

In the second phase, the Mars-centered inertial location of the lander in the interme- 
diate orbit is computed, after the landing site has been specified. The mission profile 
provides the inertial location of the orbiter in the Mars parking orbit at the time of the 
lander insertion into the intermediate orbit. In the third phase, an optimum multi- 
impulse rendezvous program, which incorporates Lawden's primer-vector theory and 
a modified Davidon search routine, is used to compute the number and location of the 
maneuvers necessary to effect rendezvous of the lander with the orbiter. This program 


optimally chooses (1) the phasing configuration for the two spacecraft before the first 
impulse, (2) the position and times for all intermediate impulses, and (3) the transfer 
time for the rendezvous. The velocity requirement for ascent (assuming inplane 
launches) is computed, and the rendezvous requirement then is plotted as a function of 
stay time on the Mars surface to produce the overall launch-window profile for the 
lander. A detailed description of the sequence of events from launch to rendezvous on 
the nominal launch day also is presented. 

As an aid to the reader, where necessary the original units of measure have been 
converted to the equivalent value in the Systeme International d' Unites (SI). The SI 
units are written first, and the original units are written parenthetically thereafter. 


SYMBOLS 


A lander cross-sectional area 

C magnitude of the intermediate impulse during lander ascent 

Cp lander drag coefficient 

G gravitational field f(x), a function of position only 

g(Z Q , t Q ) initial boundary function 

H Hamiltonian on a coasting arc 

h(Zj, t^) terminal boundary function 


h spacecraft altitude 

h apoapsis altitude 

cL 

hp periapsis altitude 

I lander ascent- engine specific impulse 

J impulsive cost, sum of the magnitudes of the applied impulse velocities 

for the N impulses 

N number of impulses 

q dynamic pressure 

r apoapsis radius of an orbital ellipse 

3 . 
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r periapsis radius of an orbital ellipse 

P 

T lander thrust 

T two-impulse Lambert reference trajectory 

t time 

V velocity 

V c characteristic velocity 

V D radial velocity 

V T tangential velocity 

W lander weight 

X displacement 

Z augmented state vector (composed of the radial and tangential velocities, 

the radial and angular position, and the corresponding Lagrange multi- 
pliers) of the lander 

z vector composed of the components of the intermediate position vector, 

the time of intermediate impulse, the initial time, and the final time 

(i pitch angle 

y vernal equinox 

At change in time 

AV change in velocity 

6( ) small deviation or perturbation of parameter ( ) 

r] true anomaly 

9 misalinement in line of apsides 

X primer vector given by Lawden's optimality condition 

v unit vector in the direction of the intermediate impulse during lander ascent 

t orbital period 

4> phase angle 
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Subscripts: 

A refers to spacecraft A 

B refers to spacecraft B 

f refers to the final condition or state of the parameter 

m refers to the condition or state at some point between the initial and final 

condition or state of the parameter 

o refers to the initial condition or state of the parameter 


Superscripts: 

’ refers to a parameter on the perturbed trajectory 

+ indicates a parameter value immediately after time t 

indicates a parameter value immediately before time t 

Operators: 

(’ ) derivative of the parameter with respect to time, d/dt 

v gradient 

REFERENCE MISSION 


A heliocentric plot of the 1981 Mars -stopover/ Venus -swingby mission used for 
the analysis presented in this report is shown in figure 1. The Earth-departure se- 
quence (computed by using the multi-orbit injection technique discussed in ref. 1) 
begins with lift-off on November 7, 1981, with the trans-Mars injection maneuver 
occurring 67 hours later. The Earth- to- Mars flight time is 266 days. The Mars orbit 
insertion (MOI) maneuver places the spacecraft in a 370- by 11 708-kilometer (200 by 
6322 nautical miles) orbit inclined 75° with respect to the Mars equator. The landing 
operations begin 20 days after MOI and result in a landing at 69° N 130° E. The 
nominal surface stay time is 30 days, and the total orbital stay time is 80 days. The 
parking orbit is designed to provide proper alinement for Mars departure (ref. 2). 

The Mars-to-Venus flight time is 127 days, with the point of closest approach to 
Venus being 4722 kilometers (2550 nautical miles). The Venus -to -Earth flight time 
is 169 days, and the Earth-entry velocity at an altitude of 121 920 meters (400 000 feet) 
is 12 048 m/sec (39 530 ft/sec). 
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180 " 



I 

Vernal 
equanox. y 


® Depart Earth Nov. 10. 1981 
© Arrive Mars Aug. 5. 1982 
© Depart Mars Oct. 24, 1982 
@ Swing by Venus Feb. 27. 1983 
(5) Arrive Earth Aug. 15. 1983 
Spacecraft (lander and orbiterl 

Figure 1. - Heliocentric view of a 1981 
80-day Mars-stopover/Venus- 
swingby mission. 


DISCUSSION OF NUMERICAL RESULTS 


Launch to Intermediate Orbit 

The nonlinear ordinary differential 
equations of motion for the time-optimum, 
constant-thrust, two-dimensional lander 
launch- trajectory computation is contained 
in reference 3. The control variable is 
the thrust orientation with respect to the 
local horizontal. The Mars-atmosphere 
model used for the analysis is the mean 
density profile presented in reference 4. 

The lander drag coefficient C D and cross- 

sectional area A were assumed to be con- 
stant. The Euler- Lagrange (costate) 
differential equations can be derived from 
equations presented in reference 3. The 
method of perturbation functions (ref. 3) 
was used to generate the numerically con- 
verged solutions. This technique, which 
is an indirect method based on the condi- 
tions required for mathematical optimality, 
is designed to satisfy these conditions by 
various iterative procedures. The bound- 
ary conditions for the state and costate 
differential equations can be expressed as 
g(Z Q ,t o ) = 0 and h(Z f ,t f ) = 0, where t Q 

and tj represent the initial and final times. 


respectively, and Z is an augmented state vector composed of the radial and tangen- 
tial velocities, the radial and angular position, and the corresponding Lagrange multi- 
pliers. The numerical integration of the state and costate equations with some Z Q 

generally will produce an h(Z^, t^) that is nonzero. The error in the terminal con- 
straint vector must be related mathematically to the initial conditions to begin the 
iteration procedure that will converge on the desired solution. The minimum-norm 
correction procedure derived and discussed in reference 5 was used to produce the 
converged solutions presented in this report. 


For the data presented, the lander ascent-engine specific impulse I was fixed 

at 380 seconds, and the initial lander ballistic parameter W /(C-.A) was chosen to be 
3 2 2 o D 

19. 15x10 N/m (400 lb/ft ). The apoapsis altitude h of the intermediate ellipse 

was always 185. 2 kilometers (100 nautical miles). 

The characteristic velocity V for the time-: optimum (free argument of peri- 

c 

apsis) transfer from the termination of a vertical rise to a 91-meter (300 feet) altitude 
to the specified intermediate orbit is presented in figure 2 as a function of the true 
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Figure 2. - Lander characteristic velocity as a function of true anomaly. 
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(c) 


h p - 121 920 meters (400 000 feet). 


91 440, 106 680, and 121 920 meters 

(300 000, 350 000, and 400 000 feet) were 

considered. Lower periapsis altitudes 

would have produced undesirable orbit- 

decay times for the lander intermediate 

orbit. It should be noted from these curves 

that the performance (measured by V ) is 

c 

insensitive to variations in true anomaly 
between ± 10°. Consequently, the choice 
of the true- anomaly angle at the 
intermediate- orbit- insertion point can be 
arbitrary within + 10°. With t ] = 5°, the 


Figure 2.- Concluded. 


characteristic velocity was plotted against 
T/W q in figure 3. The minimum V c oc- 
curs for the lowest periapsis altitude be- 
cause the apoapsis altitude is fixed at 185. 2 kilometers (100 nautical miles), and the 
minimum of that curve occurs for 0. 8 £ T/W q ^ 1.0. Based on these data, the 

periapsis altitude of the intermediate ellipse was chosen to be 91 440 meters 
(300 000 feet) (for tj = 5°), and the initial T/W q was 1. The ascent-parameter pro- 
files and dynamic-pressure profile for this case are illustrated in figure 4. 
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Figure 3. - Lander ascent performance as T " 1,e sec 

shown by the characteristic velocity as 

a function of the initial thrust-to-weight (a) Radial and tangential velocity as a 
ratio. function of time. 

Figure 4. - Typical ascent profile. 



Time, sec 

(b) Altitude and pitch angle as a function (c) Dynamic pressure as a function 
of time. of time. 

Figure 4. - Concluded. 

Geometry of an Inplane Launch 

The assumption of an inplane launch simplifies the algebra required for the 
analysis and, in most cases, produces a lower velocity requirement for the lander to 
rendezvous with the orbiter. The inplane- launch assumption provides two launch 
opportunities per day, one toward the north and one toward the south. Spherical 
trigonometry may be used to compute the Mars- centered inertial location of the lander 
at intermediate- orbit insertion once the landing site has been specified. The mission 
profile provides the inertial location of the orbiter in the Mars parking orbit at the 
time of the lander insertion into the intermediate orbit. 
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Dynamic pressure, q, Ib/ft 


Initialization of the Rendezvous Phase 


After the lander has inserted into an intermediate orbit, the rendezvous phase 
begins. The orbital orientations for the northerly and southerly launches differ, and 
the initial phase angle (central angle between the lander and the orbiter) at lander in- 
sertion depends on the day and time of launch. The nodal and apsidal rates are incor- 
porated in the orbiter equations of motion until the launch day. For a specific launch 
day, the nodal and apsidal locations of the orbiter are fixed for both daily launches 
(two-body motion); these perturbation effects are assumed to be negligible for the 
rendezvous analysis. 


Method of Solution 

The basis for the N-impulse open-time rendezvous problem originated when the 
necessary conditions for an optimum trajectory were derived by Lawden (ref. 6). This 
work was extended by Lion and Handelsman (ref. 7) to nonoptimum trajectories with the 
development of a method to determine when a reference N-impulse trajectory could be 
improved in cost by an (N + 1) impulse trajectory. Cost is defined as the sum of the 
magnitudes of the applied impulses and is discussed in the appendix. 

A gradient-type differential cost function that indicates the way to improve the 
reference solution for a fixed-time, fixed- state problem is presented in reference 7. 
These results were combined by Jezewski and Rozendaal (ref. 8) with a modified 
Davidon multivariable search routine (refs. 9 and 10). This multivariable search rou- 
tine uses a one- dimensional search (ref. 11) as well as an accelerated gradient method 
to find a local extremal. An extension of the theory to an N-impulse, fixed-time prob- 
lem that allows the optimization method to choose the true anomaly for leaving the 
initial orbit and for arriving on the final orbit is presented in reference 12. For this 
report, a further extension to the theory has been made (discussed in the appendix) 
that results in an open-time rendezvous-type solution. With a specified reference 
trajectory between an active spacecraft A and a passive spacecraft B, this method 
optimally chooses the phasing configuration of the two spacecraft, the positions and 
times of all intermediate impulses, and the transfer time needed for the rendezvous. 

A local minimum is ensured by this procedure. To guarantee that no solutions with 
extremely long transfer times would be chosen, a maximum total time boundary 

(t . . > ■ i» , ... + t, f ) of 24 hours was imposed on the solutions, 

coast in initial orbit transfer' 

It should be noted that this procedure yields only local minimums and not global 
minimums. This characteristic is inherent in the search procedure that is used, which 
assumes a unimodal function. (Unimodality is defined in ref. 13 as "in precise mathe- 
matical terms, functions y of a single variable x which, roughly speaking, have only 
one hump in the interval to be explored. ”) To overcome this difficulty, a series of 
solutions had to be obtained — each one a local minimum in a different area of the 
N-impulse solution space. 


Application of the Method 

To explain this concept further, spacecraft A and B must be considered in the 
respective coplanar ellipses with coinciding lines of apsides. If the spacecraft were 
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allowed to move in orbit to the optimum point or points for arrival and departure 
(optimal fixed-time trajectories with free true anomalies), the search procedure would 
arrive at four types of extremals — periapsis-to-periapsis transfer, periapsis-to- 
apoapsis transfer, apoapsis-to-periapsis transfer, and apoapsis-to-apoapsis transfer. 

With the same problem treated as a rendezvous problem (positions on initial and 
final orbits are linked by time) with a fixed transfer time, the cost of the solutions 
would depend on the phasing of spacecraft A and B at the beginning of the problem. If 
spacecraft A were allowed to coast At^ hours in orbit to periapsis, spacecraft B 

would have to move At^ hours plus the fixed transfer time in its orbit. The position 

of spacecraft B in orbit is not likely to be near one of the optimum places for 
rendezvous. 

Addition of the open-time capability (discussed in the appendix) to the rendezvous 
technique represents another problem. The open-time problem with an upper bound on 
total time may permit spacecraft A to coast to an orbital apsidal point at At^ hours 

and to vary the transfer time so that spacecraft B may seek an apsidal point on its or- 
bit at time At„. Thus, the transfer time for the maneuver will be At„ - At. for 
B BA 

Atg > At If the period Tg of the spacecraft B orbit is much longer than the period 

of the spacecraft A orbit, multiple extremals may be found. For example, if 

At . + r . < At„, then another solution will be found with transfer times of 
A A B’ 

Atg - (At^ + t^). These solutions form an envelope from which the minimum-cost 

trajectory between the two extremal points may be selected. The time increments 

At. and At„ are measured from the reference time. 

A B 


Nominal Results 

The nominal day of launch for the selected reference mission will now be con- 
sidered. The northerly launch will be examined first. The respective orbital orienta- 
tions of the lander and orbiter at the time of the lander orbit insertion are shown in 
figure 5. It should be noted that the lander orbit is nearly circular while the orbiter 
orbit has a 0. 5996 eccentricity. The misalinement of the line of apsides is 43. 5°. The 
near-circular lander orbit does not have well-defined periapsis and apoapsis extremals. 
This observation is supported by the results of this study, which indicate that, by using 
this method, the lander should leave at a point on its ellipse so that the total transfer 
angle between the lander (at the firing of the first impulse) and the orbiter (at rendez- 
vous) is approximately 360°. Some solutions to periapsis and apoapsis of the orbiter 
orbit were found with 180° transfers, but these solutions were similar in cost to data 
already presented and sometimes degenerated into solutions already found. In most 
cases, the generated solutions were three-impulse solutions, but some two-impulse 
solutions were found. In figure 6, the impulsive AV required for the northerly launch 

is shown as a function of the amount of time At . the lander would remain in orbit be- 

A 

fore firing the first impulse for rendezvous. The curves are numbered to differentiate 
types of solutions. Curve 1 represents the envelope of open-time rendezvous solutions 
in which the lander leaves its orbit at approximately 43° and rendezvouses with the 


$ = 43.4° imisalinement in line of apsides* 



Figure 5. - Initial orbital orientation 
for a northerly nominal launch. 


Final orbital characteristics 
h p = 370 km < 200 n. mi. J 
h a - 11 708 km (6322 n. mi.l 
Period = 7.73 hr 
Initial orbital characteristics 
hp 8 91 km (49 n. mi. i 
h a 8 185.2 km <100 n. mi. > 

Period = 1. 78 hr 

Launch A V = 4267 m/sec (14 000 ft/sec> 



x 10 


3 


o 



Figure 6. - Lander coast time in initial 
orbit before the first impulse for a 
northerly launch. 


orbiter when the orbiter first reaches periapsis 3. 04 hours after lander insertion. The 
individual points on this curve are AV minimums for the lander leaving its orbit the 
first time of arrival at the departure point, the second time, and so on. The other 
curves also are envelopes of this type. Curve 2 is the envelope of solutions in which 
the lander leaves its orbit at approximately the same point as before but does not 
rendezvous until the orbiter has reached periapsis the second time at 10. 77 hours. 
Curve 3 is the same type of envelope, except that the rendezvous occurs after the 
orbiter has reached periapsis a third time (18. 50 hours). The maximum time boundary 
of 24 hours will not permit enough time for the orbiter to reach its extremal point 
again. Curve 4 is the envelope of open- time rendezvous solutions from the lander or- 
bit to the third time the orbiter reaches its extremal at apoapsis. It should be noted 
from the transfer-to-periapsis envelopes that the minimums of all three curves have 
approximately the same AV. The difference exists only in the phasing, but, because 
the lander orbit is nearly circular, the AV penalty for phasing differences is small. 
The same effect is observed for the apoapsis curves. Thus, only one curve is 
presented. 

The initial orientation of the lander and orbiter orbits is shown in figure 7. Data 
similar to those in figure 6 are presented in figure 8 for the southerly launch on the 
nominal launch day. The identification of the envelopes is the same as before, but the 
lowest AV for the southerly launch is somewhat higher than the northerly launch 
minimum. 
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6 = 166.5° (misalinement of line of apsides) 



Periapsis of 
lander orbit 


Figure 7. - Initial orbital orientation for 
a southerly nominal launch. 


The lower of the two minimums for 
the nominal launch is chosen, and that so- 
lution is used for the nominal launch pro- 
file. This profile is presented in some 
detail. The impulsive AV needed for the 
rendezvous phase is 980 m/sec 
(3214 ft/sec). The time from insertion to 
rendezvous is 18. 5 hours. In figure 9, the 
altitudes of the lander and the orbiter are 
plotted from after the first impulse to 
rendezvous. The rendezvous- range time 
history from lander insertion to rendez- 
vous is shown in figure 10; a time history 
of the phasing from lander insertion to 
rendezvous, in figure 11; and the effect 
of each impulse on the lander orbit until 
rendezvous, in figure 12. In figure 12, 
the first impulse is shown to alter the line 
of apsides of the lander orbit to within 
2. 6 of the orbiter orbit and to raise the 
apoapsis altitude by 11 088 kilometers 
(5987 nautical miles). The second impulse 
further alines the lines of apsides and 
raises the periapsis of the lander orbit to 
370 kilometers (200 nautical miles). The 
third impulse completes the rendezvous. 
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Period -■ 1.- 78 hr 

Launch AV - 4267 m/sec (14 000 ft/sec) 
Maximum total time (coast + transfer) ■ 24 hr 
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Figure 8. - Lander coast time in initial 
orbit before the first impulse for a 
southerly launch. 



Time from rendezvous initiation, hr 


Figure 9. - Altitude time history for 
the terminal phase of the 
rendezvous. 


11 



Range, km 


1 


Time is referenced to lander orbit insertion after launch 
from Mars surface ^3 
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Figure 10.- Rendezvous- range time 
history. 
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Figure 11.- Rendezvous- phasing time 
history. 


6 * misalinement = 2.6° 

r (lander) = 3515 km (1898 n. mi. I 
P 

r (lander) ■ 14 684 km (7929 n. mi.) 
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r p (orbiter) = 3783 km (2043 n. mi.) 
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6 = misalinement = 0° 


r (lander) = 3783 km (2043 n. mi.) 

P 

r (lander) = 14 877 km (8033 n. mi.) 
a 



(a) Orbital configuration after the (b) Orbital configuration after the 

first impulse. second impulse. 

Figure 12.- Orbital- transfer geometry. 
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Impulsive AV for rendezvous, m/sec 
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Lander and 
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(c) Orbital configuration after the 
third impulse. 

Figure 12.- Concluded. 
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Data for Total Launch to Rendezvous 

The daily minimum AV for launch 
from the Mars surface to rendezvous dur- 
ing the nominal launch day has been found, 
and launches for other days of the launch 
window now will be investigated. Because 
each northerly launch has the same line of 
apsides for the lander orbit, the only 
change in the alinement for northerly 
launches for all days of the launch window 
is caused by the regression of the nodes 
and by the apsidal rate of the orbiter orbit. 
The daily alinement change for the south- 
erly launches is caused by the same effect. 
Also, each launch day (northerly and south- 
erly launches) has a different initial- 
phasing configuration than any other day. 
Thus, because the launches are inplane and 
because the lander orbit is the same for 
each day, it should not be assumed that one 
launch- window AV will be constant for all 
northerly launches and another AV will be 
constant for all southerly launches. Dis- 
cussion of the phasing for the nominal 
launch day indicates, that there should be 
variation in cost based on the initial 
changes in phasing and some based on the 
apsidal change from day to day. The cost 
(excluding launch AV) for northerly and 
southerly launches for the total Mars 
lander launch window is presented in fig- 
ure 13, which represents the results of 
this analysis. For comparison, the AV 
for the coplanar, alined line-of-apsides 
case was calculated to be 970 m/sec 
(3181 ft/sec). The southerly launch AV 
is greater than that for the northerly 
launch because of the greater misaline- 
ment of the line of apsides. 


Figure 1 3. - Rendezvous AV as a 
function of launch day. 


h 


13 


CONCLUSION 


Use of an optimum finite-thrust launch program and an optimum multi-impulse 
rendezvous program has made it possible to determine a quasi-optimum lander AV 
profile for the launch window of the 1981 Venus- swingby Mars mission. Inplane 
launches were assumed to facilitate computation, but phasing considerations and orbital 
apsidal misalinements were incorporated in the study. Results were presented in de- 
tail for the nominal launch day, and the optimum rendezvous for that day was discussed. 
The total launch-to- rendezvous AV data indicate a fairly constant rendezvous AV of 
approximately 5334 m/sec (17 500 ft/sec) with small variations caused by phasing and 
apsidal misalinement. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, June 1, 1971 
975-11-00-00-72 
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APPENDIX 


DEVELOPMENT OF OPEN-TIME RENDEZVOUS CAPABILITY 


In the calculation of optimum N-impulse trajectories for open- time rendezvous, 
a reference solution must first be established, the cost of which may be improved by 
perturbations. The impulsive cost J is defined as the sum of the magnitudes of the 
applied impulse velocities for the N impulses. 


N 

J= ^| AV i| (Al) 


The reference solution T is assumed to be the two-impulse Lambert reference 
trajectory between orbits A and B that is formed by applying impulses AV Q and AV^ 

at the times t and t., respectively. A perturbation in position 6X is applied at 


time t (t < t < t.) to the reference 

trajectory, which results in the perturbed 
trajectory T' with velocity perturbations 

6V ” and 6 V* where the minus indicates 
o 1 

an event immediately before time t and 
the plus indicates an event immediately 

after time t. The 6 V and 6 V* are 

o f 

needed to consider the variations in the 
initial and final orbits, respectively, 
caused by the perturbation. 

For the rendezvous open-time prob- 
lem, the velocities on the ends of the 
perturbed orbit are as follows (fig. A-l). 



Figure A-l. - Generalized rendezvous- 
problem geometry. 


v *" = v " + V "dt 

O 0 0 0 

(A 2a) 

V ,+ = V + + 6V +V ,+ dt 

0 O 0 O 0 

(A2b) 

V = v f~ + 5V f" + V dt f 

(A2c) 
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(A2d) 


V ' 
V f 


+ 



+ V x dt f 


The V dt terms account for the movement in the orbits and for the change in the ar- 
rival and departure times related to this movement. 

Thus, the differential cost between the reference and perturbed trajectories is 


dJ = J’ - J 


(A3) 


where 


J = 


V + - v - 

o o 


v f + - v f - 


(A4a) 


J’ = 


v * + - v * 

o o 


6V - 6V 
m 


m 


+ V 


r + 


V ' 
V f 


(A4b) 


where 6 V - 6V 

I m m 

conics of the perturbed solution at time t 


is the absolute value of the difference in velocity on the two 

By assuming small perturbations in a 


m 


gravitational field of the form G = f(x), where x is a position, and by retaining first- 
order terms, equation (A3) becomes 


dJ = X • 6V + - X. . 6V ” + <5V + - 6V “ (A5) 

o o f f m m ' ' 


where X q and X^ are the values of the primer vector X given by Lawden' s optimal- 
ity condition. Lawden' s optimality condition states that at the time of an impulse, the 
primer vector is alined with the impulse and has a constant magnitude. This constant 
magnitude is defined arbitrarily as unity by making the boundary conditions 


l (‘o) - X o - iV o/l aV ol 

(A6a) 

H*f) ■ X J - AV f/l aV £l 

(A6b) 
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The adjoint equation (a relationship linking the perturbations to the adjoint variables 
(X,-X)) is 


(X • 6V) - (X • 6X) = constant (A7) 

A method for computing X, the primer-vector derivative, is contained in reference 7. 

By use of equation (A7), evaluated at the boundaries of the two conics of T', 
equation (A5) may be transformed into 


dJ = X T • 6X - X. • 6X. + X • 6V - X 

o of f m mm 


6X 


m 


-X + -8V + +X + • 8X + + 

m m m m 


6v - 6V 

m m 


(A8) 


Because X and X are computed on the reference and, therefore, are continuous, 

X + = X " = X and X + = X = X . Also, 6X + = 6X ~ for continuity of the 
mmm mmm’mm 

position vector at t . Therefore, equation (A8) becomes 


dJ = X • 5X + - X, • 6X/ + C(1 - X • v) 
o of f v m / 


(A9) 


where C is the magnitude of the intermediate impulse and v is a unit vector in its 
direction. The arrival and departure times of the original orbits control the end-point 
motion through the equations 


dX = 6X + V dt 


(AlOa) 


or, for this discussion 


6X + = dX - V + dt (AlOb) 

o o o o 

and 

6X f “ = dX f - V f "dt f (AlOc) 
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because 


dX = V ~ dt (Alla) 

o o o 

dXf = V f + dt f (Allb) 

6X + = -(V + - V ”)dt (Allc) 

o \ o o ) o 7 

6X f “ = (v f + - V f ")dt f (Alld) 

When this equation is substituted into equation (A9), the cost differential becomes 

d J = C(1 - * m • ») - V • (V - vK - V ' (V - V)dt f (A12) 

or 

dJ = C(1 - X m • u) - |AV 0 |(* 0 + • X 0 )dt 0 - |AV f |(x ( - • Xj)dt f (A13) 

The conditions for improvement of the reference solution are obtained by inspection of 
equation (A13). 

1. If X at t Q and t^ is zero or orthogonal to and and if |X m | > 1, the 

solution will be improved by adding an impulse in the direction of X at t . 

mm 

2. If X at t and t. is not orthogonal to X and X. and is nonzero and if 

of of 

|X m | > 1, the solution can be improved by perturbing the initial and final conditions and 
by adding an impulse in the direction of at t . 

3. If X at t Q and t^ is not orthogonal to X q and X^ and is nonzero and if 
|X m | < 1, the solution can be improved by perturbing the initial and final conditions. 

Condition 2 is the most general of the conditions; conditions 1 and 3 are subclasses of 
condition 2. 
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The reference trajectory now is assumed to be a three-impulse trajectory T 
that is perturbed to obtain a T\ The differential cost (developed similarly to eq. (A9)) 
is 


dJ = X + • 6X + - X~ 
o of 


<5X r 


m 


6X 


m 



+ A 


m 


6X 


m 


6V - 6V 
m 


m 


(A14) 


For the three-impulse reference, the primer vector A m is chosen to be a unit 

vector in the direction of the intermediate impulse (6 V + - 6V ”) so that the fourth 

and sixth terms of equation (A14) vanish. The position vector X m must be continuous, 

and, because the time at which the intermediate impulse is applied varies between the 
perturbed and reference trajectories, the total change in X m is 


dX = 6X + X dt 
m m mm 


(A15a) 


dX = 6X + + X + dt 
m m mm 


(A15b) 


When equation (A15) is substituted into equation (A14), the cost function becomes 


dJ = A + • 6X + - X. - • 6X„" + (A + - A ~) • dX 
o o f f \ m m / 

- (a + *X + - a "'X ")dt 
V m m m m / m 


m 


(A16) 


By use of equation (All), equation (A16) becomes 


dJ 


= -A + • (V + - V idt - A ” • (v + - V.idt, 
o \o o/o f Vf f/f 

+ (A + - A ") • dX -(a + -X + - A " • X ") 
\m m ) m\m m m my 


m 


(A17) 
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By use of equation (A6), equation (A17) becomes 



X |AV 


dt o- 


XjAV. 


dt f + 


(* 




m 


- X 


m 


dX 


m 



X 


m 


dt 


m 


(A18) 


The Hamiltonian on a coasting arc is 


H = X • V - X • V 


(A19) 


equation (A19) now becomes 


dJ = -X 


o 


+ 



' X f| iV f| dl f + (*m + 



dX 


m 




dt 


m 


(A20) 


equation (A20) then can be written in the form 


dJ = vJ • dz 


(A21) 


where v is the gradient operator and the vector z for three impulses is composed of 
the components of the intermediate position vector, the time of intermediate impulse, 
the initial time, and the final time. Further discussion on the method of minimization 
is found in reference 10. 
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